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1. A recursive discrete Fourier transformation method 

wherein dat a values x( t) . x(t+l) , x(t+2) . x(t + 3) x(t+N-l) . 

5 x(t+N) sampled at times t, t+1. t+2. t+3 t+N-1. t+N (N 

is a positive integer which is 1 or more) each having an equal 
interval are supplied and with such N data values supplied 
since time t as a data stream, a frequency component, which 
is degree k (k is 0 or a positive integer smaller than N) obtained 
10 by carrying out complex Fourier transformation on the data 
stream, is obtained such that a real part X r (k. t) thereof 
and an imaginary part X ± (k. t) thereof are complex Fourier 
coefficients, the method comprising: 

a first step of storing the data stream x(t) , x(t+l) , 

15 x(t + 2). x(t + 3) x(t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients X r (k, t) and Xi(k. t) of the data stream stored 
temporarily at the first step; 
20 a third step of storing the complex Fourier coefficients 

X r (k, t) and Xi(k. t) obtained at the second step into the 
second memory means' temporarily; and 

a fourth step of by using data value x(t+N) supplied 
at time t+N. data value x(t) stored in the first memory means 
25 temporarily and the complex Fourier coefficients X r (k, t), 
Xi(k,t) stored in the second memory means temporarily, 
obtaining complex Fourier coefficients X r (k, t+1) and X ± (k. 
t+1) to the data stream supplied since time t+1 with respect 
to a positive constant value A for giving an amplitude value 
30 corresponding to a difference between the x(t+N) and the x(t) 
according to a following equation: 



X r {k,t + l)={x r {k,t) + ±[x{t + N)-x{tj^x 
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2. A recursive discrete Fourier transformation method 
according to claim 1 wherein inverse discrete Fourier 
transformation data subjected to inverse discrete Fourier 
transformation using a positive constant a is supplied so as 
to obtain a complex Fourier transformation coefficient to the 
supplied inverse discrete Fourier transformation data, 

the discrete Fourier transformation being carried out 
using the positive constant A corresponding to the positive 
constant a. 

3. A recursive discrete Fourier transformation method 
according to claim 1 wherein k indicating a degree of the complex 
Fourier coefficient is set to a desired value and Fourier 
transformation to the set value k is repeated at every time 
having an equal interval. 

4. A recursive discrete Fourier transformation method 
wherein with sampling frequency as f s . data values x(t ) . x(t+l) , 

x(t + 2). x(t+3) x(t+N-l). x(t+N) sampled at times t. t+1. 

t+ 2. t + 3 t+N-1. t+N (N is a positive integer which is 

1 or more) each having an equal interval are supplied and with 
such N data values supplied since time t as a data stream, 
a frequency interval given with a minimum frequency f 1 and 
a maximum frequency f 2 to the data stream is regarded as a 
measuring frequency interval, a frequency interval obtaxned 
by dividing that measuring frequency interval by the N is 
assumed to be an analysis frequency interval and a result of 
the frequency analysis obtainedby carrying out complex Fourier 
transformation at every analysis frequency interval is 
obtained as a frequency component which is k ( k is 0 or a positive 
integer smaller than N) times the frequency interval, such 
that a real part X r (k. t) thereof and an imaginary part X^k, 
t) thereof are complex Fourier coefficients, the method 
comprising: 

a first step of storing the data stream x(t) , x(t+l) . 
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x(t+2), x(t+3) x(t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients X r (k. t) and Xi(k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients 
X r (k. t) and X ± (k, t) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using data value x(t+N) supplied 
at time t+N. data value x(t) stored in the first memory means 
temporarily and the complex Fourier coefficients X r (k, t), 
X ± (k,t) stored in the second memory means temporarily, 
obtaining complex Fourier coefficients X r (k, t+1) and Xi(k. 
t+1) within the frequency interval with the minimum frequency 
f 1 and the maximum frequency f 2 to the data stream supplied 
since time t+1 with respect to a positive constant A for giving 
an amplitude corresponding to a difference between the x(t+N) 
and the x(t) according to following equations: 



X r {k,t+l) = {x r {k,t) + ±[x{t + N)-x{t)] 
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5. A recursive discrete Fourier transformation method 
25 according to claim 4 wherein inverse discrete Fourier 
transformation data subjected to inverse discrete Fourier 
transformation using a positive constant a is supplied so as 
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to obtain a complex Fourier transformation coefficient to the 
supplied inverse discrete Fourier transformation data, 

the discrete Fourier transformation being carried out 
using the positive constant A corresponding to the positive 
5 constant a. 

6. A recursive discrete Fourier transformation method 
according to claim 4 wherein k indicating a degree of the complex 
Fourier coefficient is set to a desired value and Fourier 

10 transformation to the set value k is repeated at every time 
having an equal interval . 

7. A recursive discrete Fourier transformation method 
wherein data values x(t) , x(t + l) , x(t + 2) , x(t + 3) x(t+N-l) , 

15 x(t+N) sampled at times t, t+1, t + 2, t+3, t+N-1, t+N (N 

is a positive integer which is 1 or more) each having an equal 
Ul interval are supplied and with such N data values supplied 

^ since time t as a data stream, a frequency component, which 

s is degree k ( k is 0 or a positive integer smaller than N ) obtained 

20 by carrying out complex Fourier transformation on the data 
stream, is obtained such that a real part X r (k, t) thereof 
and an imaginary part X ± (k, t) thereof are complex Fourier 
coefficients, the method comprising: 

a first step of storing the data stream x(t), x(t+l), 

25 x(t+2), x(t+3), x(t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients X r (k, t) and X ± (k, t) of the data stream stored 
temporarily at the first step; 
30 a third step of storing the complex Fourier coefficients 

X r (k, t) and Xi(k, t) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using the complex Fourier 
coefficients X r (k, t) - jX ± (k, t) stored in the second memory 
35 means temporarily, obtaining complex Fourier coefficients 
X r (k, t+1) - jXi(k, t + 1) to the data stream supplied since 
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time t+1 based on a transfer function expressed in a following 
equation . 



H(z) = A(l-z- N ). 
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where A is a positive constant for providing [x(t+N) - x(t)] 
5 with an amplitude. 

8 . A recursive discrete Fourier transformation method 
according to claim 7 wherein the positive constant value A 
for providing with an amplitude corresponding to a difference 

10 between the x(t+N) and the x(t) is capable of being set 
selectively with an inverse number of square root of 1 . N or 
1/N. 

9 . A recursive discrete Fourier transformation method 
15 according to claim 7 wherein k indicating the degree of the 

complex Fourier coefficient is set to a desired value and the 
Fourier transformation to the set value k is repeated at every 
time having an equal interval. 

20 10 . A recursive discrete Fourier transformation method 

wherein with sampling frequency as f s . data values x ( t ) . x ( t + 1 ) . 

x(t + 2), x(t+3) x(t+N-l). x(t+N) sampled at times t . t+1. 

t+2 , t+3 t+N-1. t+N (N is a positive integer which is 

1 or more) each having an equal interval are supplied and with 

25 such N data values supplied since time t as a data stream, 
a frequency interval given with a minimum frequency f 1 and 
a maximum frequency f 2 to the data stream is regarded as a 
measuring frequency interval, a frequency interval obtained 
by dividing that measuring frequency interval by the N is 

30 assumed to be an analysis frequency interval and a result of 
the frequency analysis provided by carrying out complex Fourier 
transformation at every analysis frequency interval is 
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obtained as a frequency component which is k ( k is 0 or a positive 
integer smaller than N) times the frequency interval, such 
that a real part X r (k, t) thereof and an imaginary part Xi(k, 
t) thereof are complex Fourier coefficients, the method 
comprising: 

a first step of storing the data stream x(t), x(t+l), 

x(t+2). x(t+3) x(t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients X r (k, t+1) - jX ± (k, t+1) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients 
X r (k. t+1) - jXi(k, t+1) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using the complex Fourier coefficient 
X r (k, t) - 3X ± (k, t) stored in the secondmemory means temporarily, 
obtaining complex Fourier coefficients X r (k, t+1) - jX ± (k. 
t+1) at the frequency interval given by the minimum frequency 
f 1 and the maximum frequency f 2 to the data stream supplied 
since time t+1 based on a transfer function expressed in a 
following equation: 

' c os [7ji p]- j sin \2tv p]- z' 1 ] 
l-2cos[2jrp]z _1 +z -2 J 



where A is a positive constant for providing [x(t+N) - x(t)] 
with an amplitude value, and 
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11. A recursive discrete Fourier transformation method 
according to claim 10 wherein the positive constant value A 
for providing with an amplitude corresponding to a difference 
between the x(t+N) and the x(t) is capable of being set 
selectively with an inverse number of square root of 1, N or 
1/N. 
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12. A recursive discrete Fourier transformation method 
according to claim 10 wherein k indicating the degree of the 
complex Fourier coefficient is set to a desired value and the 
5 Fourier transformation to the set value k is repeated at every 
time having an equal interval. 

13 . A recursive discrete Fourier transformation method 
for obtaining complex Fourier coefficients, wherein data 
10 values x r (t) + j Xi (t). x r (t+l) + j Xl (t + l), x r <t+2) + j Xi (t + 2). 
x r (t + 3) + dx^t+3). .... x r (t+N-l) + 3x ± (t+N-l). x r (t+N) + 

P 3Xi(t+N) sampled at times t. t + 1. t+2. t + 3 t+N-1. t+N 

S (N is a positive integer which is 1 or more) each having an 

equal interval are supplied and with such N data values supplied 
since time t as a data stream, a frequency component, which 
is degree k (k is 0 or a positive integer smaller than N) obtained 
HJ by carrying out complex Fourier transformation on the data 

H stream, is obtained such that a real part X r (k. t) thereof 

S and an imaginary part X ± (k, t) thereof are complex Fourier 

K 20 coefficients, the method comprising: 

gj a first step of storing the data stream x r (t) + jx ± (t) , 

*■* x r (t+l) + jxat + l) , x r (t + 2) + jxt(t+2) . x r (t+3) + jXi(t+3) 

x r ( t+N-1) + jXi( t+N-1) supplied since time t at time t+N-1 
into a first memory means temporarily; 
25 a second step of obtaining the complex Fourier 

coefficients X r (k, t) and X^k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients 
X r (k. t) and Xi(k. t) obtained at the second step into the 
30 second memory means temporarily; and 

a fourth step of by using the data value x r ( t+N ) + jXi( t +N ) 
supplied at time t+N, the data value x r (t) + jXi(t) stored 
in the first memory means temporarily and the complex Fourier 
coefficient X r (k, t) and X^k. t) stored in the second memory 
35 means temporarily and used recursively, obtaining complex 
Fourier coefficients X r (k. t + 1) and X ± (k, t+1) to the data 
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stream supplied since time t+1 with respect to a positive 
constant A for giving an amplitude value corresponding to the 
difference between the x r (t+N) and the x r (t) according to a 
following equation: 



5 ^ r (M+l)-|jir r M + ^r('+^)-*rM]}«»[2^] 

X, (k,t + 1) = {*,. (*,/) + ±[x, <f+N)-x, (*)]} cos 



|AT r M + ^k(^^)-^(0]]sin[2 



~N 



14 . A recursive discrete Fourier transformation method 
according to claim 13 wherein by using data value y r ( t ) + jYi(t) 
supplied at time t. data value y r (t+N) + jyiU+N) supplied 
at time t+N and complex inverse discrete Fourier coefficients 
y r (k. t) and Yi(lc. t) with the real part and the imaginary 
part obtained with respect to N data values y r (t) + jYiU) . 

y r (t+l) + JyxU+l) y r (t+N-l) + jyiU+N-l) supplied from 

time t to time t+N-1, 

a real part Y r (k. t+1) and an imaginary part Yi(k. t + 1) 
of each complex inverse discrete Fourier coefficient of N data 
values supplied since time t+1 are obtained with respect to 
a positive constant value B for giving an amplitude 
corresponding to a difference between they r (t+N) and they r (t ) 
as inverse discrete Fourier transformation data, according 
to following equations. 



25 Y,{k,t +l)-^,{k,t) + ±\y r {t +N)-y,{t)]* 
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y, (k,t + 1) - |y f M + }fo (' + * ) - y, (')]} «» [2^ 
- |y r M+ }br + yr (')]} [2^ 

and after the obtained inverse discrete Fourier transformation 
data is supplied, discrete Fourier transformation on the 
supplied inverse discrete Fourier transformation data is 
carried out, 

the recursive discrete Fourier transformation 
conducting discrete Fourier transformation by using the 
positive constant value A corresponding to the positive 
constant value B. 

15. A recursive discrete Fourier transformation method 
according to claim 13 wherein * indicating the degree of the 
complex Fourier coefficient is set to a desired value and the 
Fourier transformation to the set value k is repeated at every 
time having an equal interval. 

16 . A recursive discrete Fourier transformation method 
for obtaining complex Fourier coefficients, wherein with 
sampling frequency as fs. data values x r (t) + jXi(t) . x r (t + l) 

+ jXiCt + l), x r (t + 2) + jXiU + 2). x r (t + 3) + jXi(t + 3) 

x r (t+N-l) + jXiU+N-l). x r (t+N) + jXiU+N) sampled at times 

t, t+l. t + 2. t+3 t+N-1. t+H (N is a positive integer 

which is 1 or more) each having an equal interval are supplied 
and with such N complex data values supplied since time t as 
a data stream, a frequency interval given with a minimum 
frequency fl and a maximum frequency £2 to the data stream 
is regarded as a specified frequency interval, a frequency 
interval obtained by dividing the specified frequency interval 
by the N is assumed to be a minimum frequency interval and 
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a result of the Fourier transformation provided by carrying 
out complex Fourier transformation at every minimum frequency 
interval is obtained as a frequency component which is k (k 
is 0 or a positive integer smaller than N) times the minimum 
frequency interval, such that a real part X r (k, t) thereof 
and an imaginary part X ± (k. t) thereof are complex Fourier 
coefficients, the method comprising: 

a first step of storing the data stream x r (t) + jXi(t) . 

x r (t+l) + jXi(t + l) . Xr(t+2) + jx ± (t+2) . Xr(t+3) + jx t (t + 3) 

Xr(t+N-1) + Jx^t+N-l) . Xr(t+N) + jx ± (t + N) supplied since time 
t at time t+N-1 into a first memory means temporarily; 

a second step of obtaining the complex Fourier 
coefficients X r (k. t) and Xi (k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the complex Fourier coefficients 
X r (k, t) and Xi(k. t) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using the data value x r ( t +N ) + jxi( t+N ) 
supplied at time t+N. the data value x r (t) + jx ± (t) stored 
temporarily in the first memory means and complex Fourier 
coefficients X r (k. t) and X i( k. t) stored temporarily in the 
second memory means and used recursively, obtaining complex 
Fourier coefficient X r (k. t+1) andX^k, t+1) in the specified 
frequency interval given by the minimum frequency f 1 and the 
maximum frequency f 2 to the data stream supplied since time 
t+1 with respect to a positive constant A for giving an amplitude 
value corresponding to a difference between the x(t+N) and 
the x(t) according to a following equation: 
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17 . A recursive discrete Fourier transformation method 
according to claim 16 wherein by using data value y r ( t ) + jyi( t ) 
5 supplied at time t, data value y r (t+N) + jyiU+N) supplied 
at time t+N and complex inverse discrete Fourier coefficients 
Y r (k. t) and Y ± (k, t) with the real part and the imaginary 
part obtained with respect to N data values y r (t) + jyi(t), 
y r (t+l) + Jy±(t+1) y r (t+N-l) + jYi(t+N-l) supplied from 

10 time t to time t+N-1, 

a real part Y r (k. t+1) and an imaginary part Y ± (k, t+1) 
of each complex inverse discrete Fourier coefficient of N data 
values supplied since time t+1 are obtained with respect to 
a positive constant vaiue B for giving an amplitude 

15 corresponding to a difference between the y r ( t+N ) and the y r ( t ) 
as inverse discrete Fourier transformation data, according 
to following equations . 

jrk_ 
N 



Y r [k,t +1) = |y r M + j;\y r {t + N)-y r cos ^2 

+ |y,M+|t 1 ('+^)-y«(0][s in [ 2 f\ 



N 



20 Y, {k,t + 1) = jy, \k,t) + l[y, {t + N)-y, (*)]} cos ^2 

- |y r M + ±M'+n)- y r (')]} sin [2^ 

and after the obtained inverse discrete Fourier transformation 
data is supplied, discrete Fourier transformation on the 
supplied inverse discrete Fourier transformation data is 

25 carried out, 

the recursive discrete Fourier transformation 
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conducting discrete Fourier transformation by using the 
positive constant value A corresponding to the positive 
constant value B. 

18. A recursive discrete Fourier transformation method 
according to claim 16 wherein k indicating the degree of the 
complex Fourier coefficient is set to a desired value and the 
Fourier transformation to the set value k is repeated at every 
time having an equal interval. 



19 . A recursive inverse discrete Fourier transformation 
s method for obtaining complex inverse Fourier coefficients 

tfJ wherein complex data values x r (t) + jXi(t), x r (t + l) + jXi(t+l), 

§ x r (t + 2) + jXiU+2). x r (t+3) + jXitt + 3) x r (t+N-l) + 

ff 15 jXiCt+N-l), x r (t+N) + jxt(t+N) sampled at times t. t + 1, t + 2. 

t+3, . . . , t+N-1. t+N (N is a positive integer which is 1 or 
more) each having an equal interval are supplied and with such 
N data values supplied since time t as a data stream, an inverse 
complex Fourier transformation component, which is degree k 
i5J 20 ( k is 0 or a positive integer smaller than N) obtained by carrying 
out inverse complex Fourier transformation on the data stream, 
is obtained such that a real part X r (k. t) thereof and an 
imaginary part Xi(k. t) thereof are inverse complex Fourier 
coefficients, the method comprising: 
25 a first step of storing the data stream x r ( t ) + jx ± ( t ) , 

x r (t+l) + jXi(t+l), x r (t+2) + jXiU+2). x r (t+3) + jXi(t + 3) 

x r ( t+N-1) + jx ± ( t+N-1) supplied since time t at time t+N-1 
into a first memory means temporarily; 

a second step of obtaining the inverse complex Fourier 
30 coefficients X r (k. t) and X ± (k. t) of the data stream stored 
temporarily at the first step; 

a third step of storing the inverse complex Fourier 
coefficients X r (k, t) and X ± (k. t) obtained at the second step 
into the second memory means temporarily; and 
35 a fourth step of by using the data value x r ( t+N) +;jx 1 (t+N) 

supplied at time t+N. the data value x r (t) + jx ± (t) stored 
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in the first memory means temporarily and the inverse complex 
Fourier coefficient X r (k. t) andX^k, t) stored in the second 
memory means temporarily, obtaining inverse complex Fourier 
coefficients X r (k, t + l) and X i( k, t + l) to the data stream 
supplied since time t+l with respect to a positive constant 
B for giving an amplitude value corresponding to the difference 
between the x r (t + N) and the x r (t) according to a following 
equation: 



X r (k,t + l) = (k,t) + \[x r {t+N)-x r (*)]} cos [2^ 

10 + M + }fc + )]} sin [ 2 f ] 

X, (k,t + 1) = (*,*) + if + ")r *,(')]} «* [2^ 
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20 . A recursive inverse discrete Fourier transformation 
method according to claim 19 wherein k indicating the degree 
of the complex Fourier coefficient is set to a desired value 
and the inverse Fourier transformation to the set value k is 
repeated at every time having an equal interval - 

21 . A recursive inverse discrete Fourier transformation 
method for obtaining complex inverse Fourier coefficients 
wherein with sampling frequency as f s . data values x r (t) + 
3 Xi (t). x r (t + l) + jx^t+l), x r (t + 2) + jXi(t + 2). x r (t + 3) + 

jXi (t + 3) x r (t + N-l) + jx i( t + N-l). x r (t + N) + jXt(t + N) 

sampled at times t. t + l. t + 2. t + 3 t + N-l. t + H (N is a 

positive integer which is 1 or more) each having an equal 
interval are supplied and with such N data values supplied 
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since time t as a data stream, a frequency interval given with 
a minimum frequency f 1 and a maximum frequency f 2 to the data 
stream is regarded as a specified frequency interval, a 
frequency interval obtained by dividing that specified 
frequency interval by the N is assumed to be a minimum frequency 
interval and a result of the inverse Fourier transformation 
provided by carrying out inverse complex Fourier 
transformation at every minimum frequency interval is obtained 
as a frequency component which is k (k is 0 or a positive integer 
smaller than N) times the minimum frequency interval, in the 
form of a real part X r (k, t) and an imaginary part Xi(k, t), 
the method comprising: 

a first step of storing the data stream x r (t) + jxi(t), 

x r (t+l) + jXiU+l). x r (t+2) + j Xi (t+2), x r (t+3) + jx ± (t + 3) 

x r (t+N-l) + jXi(t+N-l), x r (t+N) + jXiCt+N) supplied since time 
t at time t+N-1 into a first memory means temporarily; 

a second step of obtaining the inverse Fourier 
coefficients X r (k. t) and Xi(k, t) of the data stream stored 
temporarily at the first step; 

a third step of storing the inverse Fourier coefficients 
X r (k. t) and X ± (k. t) obtained at the second step into the 
second memory means temporarily; and 

a fourth step of by using the data value x r (t+N) + jx ± (t+N) 
supplied at time t+N. the data value x r (t) + DXi(t) stored 
temporarily in the first memory means and preceding inverse 
Fourier coefficients X r (k. t) andX^k. t) stored temporarily 
in the second memory means, obtaining inverse Fourier 
coefficient X r (k. t + 1) andXt(k, t+1) in the specif ied frequency 
interval given by the minimum frequency f 1 and the maximum 
frequency f2 to the data stream supplied since time t+1 with 
respect to a positive constant B for giving an amplitude value 
corresponding to a difference between the x(t+N) and the x(t) 
according to following equations: 
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22 . A recursive inverse discrete Fourier transformation 
method according to claim 21 wherein k indicating the degree 
of the complex Fourier coefficient is set to a desired value 
and the inverse Fourier transformation to the set value 1c is 
repeated at every time having an equal interval. 
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